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ABSTRACT 

In this paper we advance the simple analy tic photomet ric redshift estimator 
for Type la supernovae (SNe la) proposed by IWang (120071 ). and use it to study 
simulated SN la data. We find that better than 0.5% accuracy in ^Tphot (with 
c[(^phot — Zspec)/{^ + Zspec)] < 0.005) is possible for SNe la with well sampled 
lightcurves in three observed passbands (riz) with a signal-to-noise ratio of 25 
at peak brightness, if the extinction by dust is negligible. The corresponding 
bias in Zp^ot (the mean of Zphot — -z^spec) is 5.4 x 10~^. If dust extinction is taken 
into consideration in the riz observer-frame lightcurves, the accuracy in -Zphot 
deteriorates to 4.4%, with a bias in -Zphot of 8.0 x 10^^. Adding the g band 
lightcurve improves the accuracy in Zp^iot to 2.5%, and reduces the bias in Zphot 
to —1.5 X 10~^. Our results have significant implications for the design of future 
photometric surveys of SNe la from both ground and space telescopes. Accurate 
and precise photometric redshifts boost the cosmological utility of such surveys. 

Subject headings: distance scale - methods: data analysis - supernovae: general 



1. Introduction 



The discovery of cosmic acceleration (IRiess et al. I Il998l : iPerlmutter et al. Ill999) was 



made using Type la supernov ae (SNe la) as cosmological standard candles ( PhillipsI 



1993 



Riess. Press, fc Kirshnerlll995l ). The unknown reason for the observed cosmic acceleration 
has been dubbed "dark energy" . Obtaining the spectroscopic redshifts of SNe la is the most 
costly aspect of supernova surveys. Large photometric surveys of SNe la have been planned 
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to help illuminate the nature of dark energy. Obtaining accurate and precise photomet- 
ric redshifts (photo- ^'s) will be key in maximizing the cosmological utility of such surveys 
jHuterer et aljl2004l ). 



In this paper, we advance the simple analy tic photometric redshift estimator for Type la 
supernovae (SNe la) proposed by Wang ( 2007 1^. and use it to study simulated SN la data. 
In order to derive unbiased and accurate photo-z's, it will be important to develop a variety 
of different techniques to cross-check sensitivities to various systematic uncertainties and 
correctly calculating the covariance between distance modulus and photo- 2; when fitting for 
cosmological parameters. 

We present our method in Sec. [2] and the results in Sec. [31 We discuss and summarize 
in Sec. Si 



The Method 



2.1. The analytic photo- 2; estimator 



The analytic photo-z estimator for SNe la proposed by IWangl (120071 ) is empirical, model 
independent (no templates used), and uses observables that reflect the properties of SNe la 
as calibrated standard car idles. It was dev eloped using the SN la data released by the 
Supernova Legacy Survey (lAstier et al.l 120061 ). 



This estimator uses the fluxes in griz (or riz) at the epoch of i maximum flux to make 
an effective K-correction to the i flux. The flrst estimate of redshift is given by 



■'phot 



where gj = 2.51og(/g), r/ = 2.51og(/,.), if = 2.51og(/i), and zj = 2.5\og{f^), and fg, fr, 
fi, fz are fluxes in counts, normalized to some flducial zeropoint, in griz at the epoch of i 
maximum flux. 



Next, it calibrates each SN la in its estimated rest-frame using 

M5 = 2.51og(//^7/,), 



(2) 



where f}^"^ is the i band flux at 15 days after the i flux maximum in the estimated rest-frame, 
corresponding to the epoch of At^^'^ = 15(1 + z^-^ot) days after the epoch of i flux maximum. 



^ Based on the analytic photometric redshift estimator for galaxies proposed by I Wang. Bahcall. fc Turner 



(|l998l) . 
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The final estimate for tlie pliotometric redsliift is given by 



^phot 



(3) 



1=1 



wfiere tlie data vector p = {1, gf,rf,if, Zf,i'j,i^j, Aii^}. The coefficients Cj (i=l,2,...,8) are 
found by using a training set of SNe la witli griz (or riz, for wlii cfi C2 = 0) ligf itcurves and 



measured spectroscopic redsliifts. We use tlie jackknife technique (lLuptonlll993l ) to estimate 
the mean and the covariance matrix of q (see Sec. |3]). 



2.2. Simulation of data 



SN la lightcurves for training and tes ting were simulated using the templates and various 
routines from MLCS2k2 (jjha et al.ll2007l ). We created two sets of griz lightcurves using the 
MLCS2k2 "final" templates at random redshifts in the range 0.001 < z < 1 with S/N=25 
at peak and A, the parameter that MLCS2k2 uses to describe the variation in the SN la 
lightcurves, fixed to 0. The luminosity distance fi, is calculated from the redshift assuming a 
flat cosmology with Qm = 0.3, h = 0.65 and w = —1. The particular choice of cosmology will 
not affect the photo-z estimator as long as the training and testing lightcurves are consistent 
with each other. 

The first training set consists of 100 lightcurves and was simulated with no host ex- 
tinction. The corresponding blind test set contains 1000 lightcurves similar to those in the 
training set. 

The second training set consists of 200 lightcurves, and has Ay drawn at random from 
the "default" d istribution , an ex ponential decay with argument 0.46 mag and zero for Ay < 
mag, following I Jha et al.l (120071) . The particular form of this exponential tail is derived from 
the observation (jLiralll995l ) that SNe la have a common color at ~ 40 days past peak. We 
assume the extinction law to be Ry = 3.1. We used a larger training set of 200 simulated 
lightcurves for training to compensate for the larger variation in the data caused by including 
extinction. The corresponding blind test set contains 1000 lightcurves similar to those in the 
training set. 

The actual training sets and blind test sets used are slightly smaller than described 
above, because only SNe la with riz lightcurves, i.e. z < 0.95, are used in this paper. We 
use the training sets to derive the photo- 2; coefficients and their covariance matrix (see Eq. [3] 
and Sec. [3]), and the blind test sets to evaluate the performance of the photo- 2; estimator. 
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3. Results 

We use a slightly modified version of the jackknife technique ( LuptonlllQQsl ) to estimate 
the mean and the covariance matrix of q. From the training set containing N SNe la, we 
extract N subsamples each containing — 1 SNe la by omitting one SN la. The coefficients 
c[^^ (i=l,2,...,8) for the s-th subsample are found by a maximum likelihood analysis matching 
the predictions of Eq. [3] with the spectroscopic redshifts. 

The mean of the coefficients Cj (i=l,2,...,8) are given by 



1 ^ 

(-^) = ^E-S^- (4) 



s=i 

Note that this is related to the usual "bias-corrected jackknife estimate" for Cj, c/, as follows: 

ci = cf + (iV - 1) (cf - (q)) , (5) 

where cf are estimated from the entire training set (with N SNe la). We found that for 
small training sets (with < 20) that include SNe la at 2; ~ 0, c/ give biased estimates of 
Ci by giving too much weight to the SN la with the smallest redshift. For training sets not 
including nearby SNe la, (q) and cf are approximately equal. We have chosen to use (q) 
from Eq. [Has the mean estimates for q to avoid biased -Zphot for SNe la at z close to zero. 

The covariance matrix of Cj (i=l,2,...,8) are given by 

Cov(c., c,) = ^ Yl ^ - - (6) 

s=l 

Since Zphot = Yli^iVii ^^^phot = YliVi^^^ii since the uncertainty in Zphot is dominated by the 
uncertainty in q. Therefore estimated error in Zphot is 



8 



1/2 



d^phot = < 2^ iJPi] Cov(q, Cj) [pj] [> , (7) 
I i=i j=i 

where Pi and pj (i,j=l,2,...8) are the i-th and j-th components of the data vector p = 
{'^,9f,rf,if,Zf,ij,ij,Aii5}. 

In applying the photo- 2; estimator to simulated data, we find that applying it to the 
entire data set leads to substructure in the estimated ^phot- Therefore, we divide the train- 
ing set into three i maximum fiux ranges (corresponding to three roughly equal redshift 
intervals), and applying the photo- ;z estimator to each fiux range. 
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Fig. [T] shows the resultant performance of the photo-z estimator for simulated SN la 
lightcurves in riz, with S/N=25 at peak brightness, and zero extinction due to dust. The 
top panel shows the training set of 97 SNe la, with different point types and color denoting 
different flux ranges. The q's derived from the training set are used to predict Zphot for a 
blind-test set of 940 SNe la (divided into the same flux ranges as the training set), with the 
results shown in the bottom panel of Fig. [H 

The results for the blind-test set closely mimic that of the training set. Since the data 
are very sparse at 2; < 0.02 due to the small size of the training set, the estimated uncertainty 
on 2phot, d2;phot, becomes large for nearby SNe la. Since dzphot is a reliable indicator of the 
accuracy of -Zphot, we can use dzphot to exclude the SNe la with estimated 2;phot not suitable 
for inclusion in the cosmological data set. 

We find that dzphot < 0.1 only excludes 9 out of the 940 SNe la in the blind-test set, all 
at z < 0.02. For the culled blind-test set of 931 SNe la (d2;phot < 0.1), o"[(2phot — 2spec)/(l + 
^spec)] = 0.0043. The corresponding bias in 2;phot (the mean of -Zphot — ^igpec) is 5.4 x 10""^ overall 
and is [0.48, 0.79, 1.5, -0.81] x lO^^ ^^e 4 redshift bins bounded by [0, 0.25, 0.5, 0.75, 1]. 

Fig. [2] shows the performance of the photo-z estimator for simulated SN la lightcurves 
in riz, with S/N=25 at peak brightness, and extinction due to dust parameterized by Ay- 
Again, we cull the blind-test set by requiring that dzphot < 0.1; this excludes 16 out of the 
956 SNe la, most of which are highly extinguished or have very low z. The inclusion of dust 
extinction leads to a deteriorated accuracy in Zphot of 4.4%. The corresponding bias in Zp^ot 
is 8.0 X 10^^ overall and is [—9.3, —22., 25., 53.] x 10"^ in the 4 redshift bins bounded by 
[0,0.25,0.5,0.75,1]. 

Fig. [3] is similar to Fig. [21 but only the subsets of data (both the training set and the 
blind-test set) with g band lightcurves are included. Note that the i maximum flux ranges 
in Fig. [3] are similar to Fig. [21 but slightly adjusted to make sure that the sub sample with 
the smallest fluxes contains enough SNe la. The SNe la at 2; > 0.54 in our simulated data 
set do not have g band lightcurves, as these require rest-frame lightcurve templates bluer 
than u band, which are not yet readily available. 

Imposing d^phot < 0.1 excludes 12 out of the 537 SNe la, most of which are highly 
extinguished or have very low z. Adding the g band lightcurves improve the accuracy in 
^^phot to 2.5%. The corresponding bias in Zphot is reduced to —1.5 x 10~^ overall and is 
[-6.0,3.2] X 10-3 in the 2 redshift bins bounded by [0,0.266,0.532]. 
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4. Discussion and Summary 



Accurate and precise pho to-z's will be critical for enabling cosmology with future large 
SNe la photometric surveys (iHuterer et al.l 120041 ) . In this paper, we have advanced the 
empirical and analytic photo- 2; estimator proposed by IWangl (120071 ) and used it to study 
simulated data in order to derive the survey requirements for accurate and precise photo-z's 
for SNe la from large photometric surveys. Future work will explore the critically-important 
question of the degree of covariance between the redshift determined using this method and 
the distance modulus determined by using a luminosity distance fitter since it is these two 
parameters that form the basis for the cosmological measurements with SNe la. 

We apply the photo-z estimator to riz data divided into i maximum flux ranges that 
roughly correspond to equal redshift intervals (see Fig. [T]and Fig. [2]). We use dzphot < 0.1 
as a quality measurement to cull the blind-test set. This cut only excludes < 2% of the 
SNe la, most of which are highly extinguished or are at very low z. To derive accurate 
cosmological constraints, it is optimal to combine the very large set of SNe la with well 
estimated photo- 2;'s with a small spectroscopi c sample of SNe la at very low z expected from 
ongoing nearby surveys (IHuterer et al.ll2004l ). SNe la a.t z < 0.025 are currently excluded 
from SN la cosmological analyses even for spectroscopic data, because peculiar velocities 
due to cosmic large scale structu re modifies the Hubble diagram of nearby standard candles 
( Wang. Spergel. fc Turner 1 199 sl), and have th e large st impact on th e near es t SNe la (see, 



Zehavi et al 



(l2007h ). 



Jl998h : ICoorav fc Caldwelll J2006h : IHuI fc Green"3 J2006h : IConlev et al 



We find that determining ^phot's with better than 0.5% accuracy (RMS[Azphot/(l+^spec)] 
0.0043%) and negligible bias {{z 



phot 



■'spec 



) = 5.4 X 10~ ) is possible if dust extinction 
can be neglected, for SNe la with well sampled lightcurves in riz and S/N=25 at peak 
brightness. The inclusion of dust extinction leads to a drastic deterioration of the accuracy 
of 2;phot (RMS[Azphot/ (1 + ^spec)] = 4.4%, and (^phot -^^spec) = 8.0 x 10"^). Adding the g band 
lightcurve significantly improves the accuracy and precision of Zp^ot with RMS[A2phot/(l + 
^^spec)] = 2.5%, and (2;phot - ^^spec) = -1-5 x 10"^. 



Future su pernova surveys can easily obta i n the multi-band ph otometry of a huge number 
of supernovae (IWangll2000l : IWang et al.ll2004l : iPhillips et al.ll2006l ). We only considered griz 
photometry in this paper, due to the limitation of currently available SN la lightcurve tem- 
plates. The addition of photometry at longer wavelengths should provide more dramatic im- 
provement on the accuracy and precision of 2;phot sinc e extinction decreases with wavelength , 
and SNe la are better standard candles in the near IR (IKrisciunas. Phillips, fc Suntzeffll2004l ) . 
In future work, we will study simulated SN la data over the entire wavelength range of rel- 
evance to future observations, and examine the impact of realistic survey and observational 
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parameters on the precision and accuracy of our photo-z estimator. 

Our results are encouraging for the prospect of using photometric surveys to constrain 
cosmology, such as those planned for the Advanced Liquid-mirror Probe for Astrophysics, 
Cosmology and Asteroids (ALPACA)@; Pan-STARR^I; the Dark Energy Survey (DEsfl 
and the Large Synoptic Survey Telescope (LSSTfl. 
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Fig. 1. — The analytic photo-z estimator applied to simulated SN la lightcurves in riz, with 
S/N=25 at peak brightness, and zero extinction due to dust. 
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Fig. 2. — The analytic photo-z estimator applied to simulated SN la lightcurves in riz^ with 
S/N=25 at peak brightness, and dust extinction parameterized by Ay. Dust extinction leads 
to deterioration in the accuracy of Zphot- The y-axis range is the same as in Fig. [3l 
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Fig. 3. — The analytic photo- 2; estimator applied to simulated SN la lightcurves in griz, 
with S/N=25 at peak brightness, and dust extinction parameterized by Ay- Adding the g 
band lightcurves improves the accuracy of Zphot- The y-axis range is the same as in Fig. O 



